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Abstract 

We consider a possibility to use the solar neutrinos for studies of small 
scale structures of the Earth and for geological research. Effects of thin layers 
of matter with density contrast on oscillations of Beryllium neutrinos inside 
the Earth are studied. We find that change of the 7 Be neutrino flux can reach 
0.1 % for layers with density of oil and size 20 km. Problems of detection are 
discussed. Hypothetical method would consist of measuring the 7 Be— flux 
by, e.g., large deep underwater detector— submarine which could change its 
location. 
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1 Introduction 



In geophysics the matter density profile of the Earth is determined by mea- 
surements of the seismic wave velocity profile. Seismic method is used widely 
also for the oil search, in particular, in the decision process before a trial well 
is drilled. This method has, however, shortcomings and alternatives are 
highly appreciated. 

Enormous penetration ability of neutrinos is a source of strong tempta- 
tion to use them in tomography and geological studies [I] , [2] , , |3] The 
first proposal was based on inelastic scattering of high energy neutrinos pro- 
duced by accelerators [T]. The authors have suggested several methods which 
use difference of the neutrino interactions in matter with different chemical 
composition and density. 

Qualitatively different proposal is based on the elastic forward scattering 
(refraction) of neutrinos [2] , [3] , |I] . Refraction modifies properties of neutrino 
oscillations changing both the oscillation length and depth The effect 
depends on density of medium. Consequently, an appearance of layers with 
different densities on the way of neutrinos changes the oscillation pattern. In 
this connection, some possibilities to use the neutrino superbeams as well as 
the beams from neutrino factories have been considered [BJ. 

There is a number of studies of the Earth matter effect on oscillations 
of solar neutrinos. It was marked that the solar neutrinos can give informa- 
tion about large scale density distribution inside the Earth. In particular, 
properties of the mantle and the core of the Earth can be studied [7j. 

The Earth matter effect has been considered for the 7 Be neutrinos jH] 
in connection to BOREXINO jTUj and KamLAND ^T] experiments. It was 
found that the effect is negligible for the LMA solution, at least for these 
detectors. 

Main conditions to use neutrinos for the geological studies and searches 
for oil and minerals are: 

1) sensitivity to small scale structures: d < 100 km; 

2) possibility to move both a source and detector of neutrinos, so that 
substantial part of outer layers of the Earth mantle can be scanned. 

No realistic proposals which satisfy these conditions have been published 
so far. In this paper we will consider a possibility which is based on: 

• Detection of the solar 7 Be— neutrinos: For the LMA solution of the 
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solar neutrino problem the oscillation length of neutrinos with E^ e ~ 
0.86 MeV is in the required range L,/2 ~ (10 -j- 20) km. 

• Measuring the time dependence of the neutrino flux by deep underwater 
detector - submarine. 

The paper is organized as follows. In section 2 we consider the effect of 
thin layers of matter on neutrino oscillations (see for some related work). 
We study dependence of the effect on neutrino parameters and properties of 
layers. In section 3 we consider a possibility to identify the effect of cavities. 
In section 4 we discuss problems of detections. Concluding remarks are given 
in section 5. 

2 Effects of thin layers on neutrino oscilla- 
tions 

Let us consider a system of two mixed neutrinos with values of oscillation 
parameters Am 2 and sin 2 29 from the LMA solution region [T^j. (We will 
comment on the effect of the third neutrino at the end of this section). In 
the case of the LMA solution, the electron neutrino produced in the center of 
the sun is adiabatically converted to a combination of the mass eigenstates 
v\ and v-i which is determined by the mixing angle, (P ml in the production 
point: 

v e — > cos 6^ v\ +sin#^ u 2 ■ (1) 
The angle 9^ is given by: 

tan 2 2^ = tan 2 2fl(l- 2V J E \ ' , (2) 
V Am 2 cos29J ' v ' 

where E is the neutrino energy, V® is the matter potential for neutrinos in 
the production point. V e = \Z2GfpNaY £ , Gp is the Fermi coupling constant, 
p is the density of medium, iV^ is the Avogadro number and Y e is the number 
of electrons per nucleon; V e ° = V e (p°, F e °). 

Conversion effect should be averaged over the neutrino production region, 
and in what follows we will describe this averaging by the effective mixing 
angle 6^. 
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Due to loss of coherence j!4j . neutrinos arrive at the surface of the Earth 
as incoherent fluxes of v\ and v 2 with relative admixtures given by cos 2 (P m 
and sin 2 (P m correspondingly. 

Let us consider oscillations inside the Earth. The probability to find v e 
in the detector can be written as: 

P = cos 2 6° m P le + sin 2 6° m P 2e = cos26° m P le + sin 2 6° mJ (3) 

where P\ e and P 2e are the probabilities of v\ — > v e and v 2 — > v e transitions 
in matter of the Earth correspondingly. 

Suppose neutrinos cross the cavity with length d and density p^ which 
differs from the average density of surrounding matter, p. That is, the neu- 
trinos cross consecutively three layers: with lengths and densities (h, p), (d, 
Pd), (h, p)- The total length of the trajectory is L = l\ + d + l 2 . 

Let us find an effect of the cavity on oscillations. Introducing, P® e , the 
v \ — * v e oscillation probability in absence of cavity (it corresponds to oscil- 
lations in the unique layer of density p and length L) we can write P as: 

P = P° + AP, (4) 

P° = sm 2 e° m + cos 26° m Pl (5) 

AP =cos2^(P le -P 1 ° e ). (6) 

Here P° is the probability to find v e in the absence of cavity and AP is the 
change of the probability due to the cavity effect. 

We will calculate the effect for the beryllium neutrino line with energy 
Ese = 0.86 MeV. For such an energy the matter effect on mixing and oscil- 
lations is very weak. It is determined by a small parameter: 

e = ?M « 26- 10- ( P ) f 7 - 10 ' 5 ^ 2 ) ( M (7) 
Am 2 ' V2.7g/cm 3 ; V A ™ 2 A - 5 ^ 

which characterizes deviations of the mixing angle and the oscillation length 
in medium from their vacuum values. We will use e as an expansion pa- 
rameter and find effects in the lowest order in e. The oscillation length is 
determined by Am 2 : 

C « « 30.7 km ( 7 '™2f ) ■ (8) 
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The probability of v\ — > v e oscillations equals 

P u = |cos# S n + sinO S 12 \ 2 , (9) 

where S = | \Sij\ | is the evolution matrix in the mass eigenstate basis, v. The 
matrix S can be written as: 

S = UD 2 U ] U d D d U\UD 1 U^ . (10) 

Here U = U(p) and Ud = U (pd) are the mixing matrices of mass eigenstates 
in the layers 1, 2 and in the cavity correspondingly. They are defined in such 
a way that v = Uu m , where u m is the basis of eigenstates of the Hamiltonean 
in matter. In (jl(Jj) . the diagonal matrices, 

Di = diag(l, e m ), i = 1,2, d, (11) 

describe evolution of the neutrino eigenstates, u m , in layer i. is the phase 
of oscillations: 

/\rn /■ / 

$i = -^y/ (cos26 + sin 2 20 , ( 12 ) 
6j = e(pi). In absence of cavity we have 

S = S = UDU + , (13) 

where D = diag(l, e**) and 

Am 2 L 



$ = — ^/(cos20-e) 2 + sin 2 2# (14) 

is the total oscillation phase in absence of the cavity. 

Using Eqs. (|9I1UI11|) . the straightforward calculation leads to the V\ — > u e 
probability in the lowest order in e: 

P u « Pu + ^ sin 2 26 sin ($2 + y ) sin , (15) 

where 

r _ ( Y ep) ~ ( Y edPd) (16) 

^~ YeP 1 ' 
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is the density contrast. In absence of cavity the probability equals 



po = cos 2 9 - e sin 2 29 sin 2 - . (17) 
And according to eq. (jHJ) and ([To]) the effect of cavity can be written as 

AP = e«e cos 29° m sin 2 29 sin ($ 2 + ^) sin -± . (18) 
The following remarks are in order: 

According to eq. (fT%|) the effect of cavity is proportional to expansion 
parameter e, the density contrast, £, and cos 2^. For maximal mixing in 
the production region the effect would be zero. For the best fit point of the 
LMA solution we get cos 29^ « 0.3. 

Important feature of the result (fTSj) is that AP does not depend on the 
oscillation effect in the first layer (before cavity), in particular, it does not 
depend on the phase $i. In contrast, AP does depend on $ 2 - the oscil- 
lation phase in the layer between the cavity and detector. This property 
appears in the lowest order in e and is related to the fact that initial state 
is the mass eigenstate and final state (at the detector) is the flavor state. 
It can be shown that if initial state is v e and final state - V\, the proba- 
bility AP does not depend on $ 2 (phase between cavity and detector) but 
will depend on $i— the phase acquired in the first layer. Thus, the inter- 
change of v\ and v e in the initial and final states leads to change $2 ~ ► 3>i- 
This feature is important for identification of the cavity effect (see section 3). 

Let us estimate the size of possible effect depending on parameters of 
neutrinos and the cavity. The cavity produces the change of the event rate 
during the time when the neutrino beam crosses the cavity. The relative 
value of the change AN/N = AP/ P equals in the lowest order in e 

A/> /(AmV) efsin($ 2 + %W% ; (19) 



po J y ' ' * V 2 J 2 
where 

= 2cos2esm^ 

JK ' 1 + cos 20°, cos 20 V ' 

Here we have taken P° = (1 + cos 29® m cos 29)/ 2 neglecting the Earth matter 
effect. The dependence of f(Am 2 ,9) on tan 2 9 for Am 2 from the allowed 
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LMA region is shown in Fig. 1. According to Fig. 1 the function /(Am 2 , 9) 
has maximum f max ~ 0.4 at about the best fit value of tan 2 9. 




The relative change of the probability can be rewritten as 

AP / <I>A $^ 

A max sin $ 2 + sin , (21) 



where 

AP 

A max = — - = eU(&rn 2 ,6) (22) 

-*e max 

is the maximal value of AP/P° which corresponds to both sines being 1. The 
effect increases with decrease of Am 2 . For minimal allowed value Am 2 = 
5 ■ 10" 5 eV 2 , and tan 2 9 = 0.3, £ = 0.64 we get 

A max = 0.1%. (23) 

In Fig. 2 we show A max as the function of Am 2 and tan 2 9. The effect weakly 
depends on 9 in the allowed LMA region. 
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Fig. 2. Maximum value of the relative effect, A max = (AP/P) max , as the 
function of Am 2 and tan 2 9. 



If detector is placed underwater, neutrinos will cross an additional layer 
of water, lw, which can be rather large. Straightforward calculations lead to 
result, similar to eq. ()19)) in which 

. , Am 2 l w 

$ 2 _> $ 2 = $ 2 + $ w , $ w ~ 2 ^ , 

where $ w is the phase acquired by neutrinos in water. 

Let us comment on possible effect of the third neutrino. We consider the 
scheme in which v% is separated from two others by the mass gap Am| 3 = 
Am 2 atm = (2 — 3) • 10~ 3 eV 2 . (So, the scheme explains the atmospheric neu- 
trino data.) We assume that this neutrino has small admixture of the electron 
neutrino described by U e ^ < 0.15. The influence of matter on this mixing 
is determined by e atm = e (^ m atm) ~ ? • 1CT 5 . Consequently, variations of 
the flux due to cavity effect should be smaller than 7 ■ 10 -5 . In fact, the 
observable effects are further suppressed by smallness of U e s and averaging 
of oscillations associated with the third neutrino. For Am 2 atm the oscillation 
length is smaller than 1 km. Also interference between the modes of oscil- 
lations driven by Am 2 atm and solar Am 2 produces negligible effect. So, we 
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conclude that the cavity effect is mainly due to oscillations driven by solar 
Am 2 . 

3 Effect of cavity. Exposure time 

Let us consider the time dependence of the electron neutrino flux in, e.g., 
deep underwater detector. For fixed position of the detector, the time depen- 
dence is due to rotation of the Earth. The rotation leads to certain change 
of the trajectory of solar neutrinos inside the Earth, and consequently, to 
modification of the oscillation effect. In general, the change of signal with 
time is both due to the cavity effect and the effect of regeneration in P° (see 
(jHJ)). The problem is to disentangle these two effects. 

Let t exp be the exposure time - the time during which the trajectory of 
solar neutrinos crosses the cavity. The effect of cavity consists of an additional 
change of the neutrino flux during exposure time t exp on the top of usual 
regeneration effect. 

According to (j5J) and (|17|) the probability in the absence of cavity equals 



In principle, this probability can be predicted with high enough accuracy. 
Indeed, the average radial density profile (THj is a subject to local variations 
at the level (5 — 7)% [IB]. Therefore one expects ~ (5 — 7)% variations of 
e(t). This uncertainty can be important for very small values of the den- 
sity contrast: £ ~ (5 — 7)% or very small phases in Eq. (fT9"j) . Thus, one 
should look for the deviation of the time dependence of signal from the sig- 
nal expected according to ()24|) . Identification of the cavity effect depends 
on specific conditions of observation. If the cavity is close enough to the de- 
tector, the distance l±, and therefore L, can change significantly during the 
exposure time, so that the averaging over the phase $ will occur. As a result 



practically does not depend on time, and whole time dependence is due to 
presence of cavity. 



P°{t) = sin 2 9° m + cos 29° m cos 2 9 - e(t) sin 2 29 sin 



2 $00 



(24) 



2 




(25) 
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Due to smallness of relative change of the probability, AP/ P, one needs 
to collect the integral effect during long time of observation. Therefore, the 
effect of cavity is determined by the relative change of probability averaged 
over the exposure time: 



-^—A max f dt sin { $' 2 (t) + 



AP _ 1 „ , *M sin Ml. (26) 



In general, the cavity has complicated form, so that d = d(t) can change with 
time during exposure interval significantly. Averaging over d gives 



AP 1 -A max 



P t 

1 °exp 



J dt cos <$>' 2 (t) . (27) 



The size of the cavity should not be much smaller than the oscillation 
length. For the Am 2 ~ 7- 10~ 5 eV 2 the oscillation length is 30 km. Therefore 
maximum effect should be seen for the cavity size l u /2 ~ 15 km. For smaller 
cavities the effect will be suppressed according to the Eq. (fT§]). 

AP is the periodic function of the phase $ 2 . Averaging over $ 2 leads to 
disappearance of the cavity effect AP = 0. Suppose averaging is absent. By 
moving the detector one can vary distances l 2 and lw Selecting the distance 
li and lw in such a way that <3> 2 + $d/2 = tt/2 + irk (k = integer) one can 
maximize the effect 

AP = e f sin 2 26 cos 26° m sin -± . (28) 

Correspondingly, there are zones of 1% with high and low sensitivity. Moving 
the detector by l v /2 one can cover whole space. 

The width of the beryllium neutrino line, APe e rs 1.5 keV [T7] . is very 
small. Consequently, no averaging over the energy occurs. Averaging would 
require the length of trajectory: 1% > l^Ese/ 1 AEse ~ 2 ■ 10 4 km which is 
larger than the Earth diameter. 

The sign of the effect depends on the sign of density contrast, the distance 
between the cavity and the detector, l 2 , and the size of the cavity itself. In 
particular, if the cavity is close enough to the detector, so that $ 2 tt, and 
therefore AP ~ £sin 2 $d/2, the v e flux (signal) increases, when neutrinos 
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cross the cavity for the positive contrast £ > pd < P2 (which is satisfied 
for oil or water), and it decreases for the negative contrast. Similar situation 
is for $2 ~ 2irk, k = 1,2, .... The sign of the effect changes if the detector 
approaches or removes from the cavity by the distance ~ /„/2. 

Detecting the cavity from the first position of the detector one determines 
the line along which the cavity is situated. Observing the cavity from two 
different positions of detector one can reconstruct location of the cavity unless 
another cavity is detected from the second position. The latter can be checked 
by observation from the third position of the detector. 

In general, one should take into account finite size of the region in which 
beryllium neutrinos are produced inside the sun. The diameter of the pro- 
duction region (~ O.2i? ) has the angular size at the Earth 9 Be = 9 • 1CT 4 
radian. Therefore, the cavity with transverse size h shields the production 
region completely if 

k < A = 1100 km . (29) 

If the cross-section of the cavity, S d , is smaller than (hOBe) 2 an additional 
suppression factor appears in the size of the relative change of the flux: 

-p* (5^P (30) 

The integral effect is determined by the exposure time, which in turn 
depends on (i) the size and shape of the cavity, (ii) its location, (iii) position 
of detector and (iv) season of the year. For simplicity we will consider the 
cavity of the parallelepiped form situated in the horizontal layer, with h being 
the thickness in the vertical plane, l w — the width in the horizontal plane, and 
as before, d— the size of the cavity in horizontal plane along the direction 
from the Sun to the detector. 

Let tf^v be the exposure time during the day, that is, during one crossing 
of the cavity by the neutrino trajectory. If v— trajectory crosses the cavity 
N days during a year the total exposure time during the year equals: 

t exp ~ tt y p ■ N. (31) 

Let us estimate the exposure time for different possible positions of the 
cavity. 
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Suppose the cavity is in the equatorial region. If the detector is situ- 
ated in the same horizontal plane as the cavity, the diurnal exposure time is 
determined by the cavity thickness: 



where t day = 24 hours. If the detector is out of the plane or cavity is inclined, 
the exposure time can be larger. However, at the same time the length of v 
trajectory in the cavity will be shorter. Both factors can compensate each 
other in the integral number events. 

The cavity will shield the detector during N = 365 days, so that the total 
exposure time during the year can be as large as ~ 14 hours. 

Let us consider the cavity near the North pole. Here the search will be 
possible in spring or in autumn when the sun appears at horizon. At this 
time the sun rises by ~ 7- 10~ 3 radian/day. Therefore the sun will be shielded 
by the cavity during N — h/(7 • 10 _3 / 2 ) days. For thickness h — 1 km we get 
N = 1.5 days. The diurnal exposure is determined by the width l w , so that 
total exposure time equals 

LtdayN (lwidth\ /100km \ / h \ . . 

^ = ~2^ day " 1 h ° Ur (20km- ) [— ) { Tkm" J ' (33) 

Let us consider now the search at latitudes 67° which correspond to the 
polar circle. In this case in December the sun is at horizontal (south) direc- 
tion for 18 days a year. The diurnal exposure time is the same as in previous 
case. As a result, t exp = 13.5 hours. Similar consideration holds for observa- 
tion in June at the south latitude 67°. 



4 Detection 

Clearly, detection of so small ( < 0.1%) effect at low energies and during 
restricted time interval is a very challenging task. It may require new tech- 
nological developments. Here we will give simple estimations which help to 
understand problems of the method. 

Let us consider a possibility to detects the effect of the cavity by large 
underwater detector-submarine. 
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The background conditions should be as good as in the underground 
detectors of Beryllium neutrinos (TU1 HI] . To suppress the background gen- 
erated by the cosmic rays, the detector should be placed deep underwater 
— at about 3500 m below sea level. This certainly restricts applications of 
the method and makes it more difficult technologically. Indeed, it is unlikely 
that, e.g., oil reservoirs are very deep although such possibility is not ex- 
cluded [18] . On the other hand, there is not too many places where the sea 
becomes very deep near the cost. 

The background and therefore the depth could be reduced for detectors 
which use light elements (e.g., helium) and efficient anti-coincidence shield- 
ing. 

Let us estimate the change of number of events due to presence of cavity 
and required size of the detector. For definiteness we will consider detection 
of neutrinos via the ve~ scattering in scintillator. The rate of events due to 
the Beryllium neutrino flux without oscillations is [TU] : 

Then in presence of neutrino conversion the expecting rate is 

r , (35) 



P + \(l-P) 



where P is the survival probability of the electron neutrinos. Here we have 
taken into account that u e e~ cross section is about 5 times larger than z/^e - 
cross section for 1 Be neutrino scattering with energy transfer (0.25 - 0.664) 
MeV. Taking typical value P° ~ 0.6 we get 

events . . 

r = 14- . 36 

kt ■ hour 

Total number of events detected during the year equals 

N = r-M d . tf:;, (37) 

where M d is the mass of detector. 

The effect of cavity consists of the increase (or decrease) of the number 
of events by amount 



A ^ = JV -^-TTT?4p5 = a7JV -^ (38) 
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(see Eq. ()26|l). Last factor in this equation reflects dumping effect due 
to contribution of neutral current interactions. This factor is absent for 
detection via the charged current interactions. 

Let us find the number of events needed to establish an excess (or deficit) 
of the number of event due to presence of the cavity at ka level (k is inte- 
ger). We take an ideal situation of 100 % efficiency of detection, absence of 
background, absence of averaging, etc.. In this case the possibility to identify 
the effect of cavity will be determined by statistical fluctuation of the neu- 
trino signal during the exposure time: AN sta t = VN- The effect of cavity 
should be larger than the fluctuation. The effect of cavity at ka level will be 
achieved if AiV c > kAN stat = kyN. From this condition we find using Eq. 
(HH): 

iV>2.A; 2 .^ 2 . (39) 

For k = 3 (3a— level) and AP/P = 0.001 this inequality gives 

N > 1.8 ■ 10 7 events. (40) 
Consequently, the mass of detector should be 

n 2-k 2 /afv 2 
Md ^ — = — ■■hH • ( 41 ) 

' v exp 1 v exp \ 1 J 

Taking the exposure time t exp ~ 10 hours we find M = 130 Mt. In general, 

. » k 2 lOhours . „„. 

M d = 130Mt — — . (42) 

" ^exp 

For the density of scintillator 0.8 g/cc a detector could have dimensions 
~(500 m) 3 . (Compare with SuperKamiokande which has diameter 39 m 
and the height 42 m.) Notice that for such a large detector the problem can 
appear with fiducial volume. In our estimations of the exposure time we have 
neglected the finite sizes of the detector. In fact, only part of the so large 
detector can be shielded by the cavity: that is, the trajectories of neutrinos 
which cross the cavity will cover only part of the detector during all exposure 
time. Therefore the effect of cavity will be further suppressed. 

We have found that the size of detector (even in the most favorable sit- 
uation) should be an order of magnitude larger than the size of future un- 
derground detectors which are now under consideration Solution of this 
problem may require next step of technological developments. 
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5 Discussion and conclusions 



We have studied effects of the thin layers of matter on oscillations of solar 
neutrinos, in particular, on the Be-neutrinos. The effect is proportional to 
e, the effective mixing parameter cos 26*^ in the production point, and the 
density contrast. 

The oscillations occur in the vacuum dominating regime and the matter 
effect appears as the small correction to the vacuum oscillations pattern. For 
the beryllium line the oscillation length is about 20 - 40 km and the matter 
effect parameter e ~ 3 ■ 10~ 3 . This parameter put the absolute upper bound 
on the effect for a given energy. The only possibility to further enhance the 
effect is to have several layers with configurations which add up their effects 
coherently |2"Uj . 

The interesting feature of the oscillation result is that for solar neutrinos 
it does not depend on the oscillation phase acquired by the neutrinos before 
the layer with density contrast. 

We have considered a possibility to study small scale structures inside 
the Earth via oscillations of the solar Beryllium neutrinos. Such a study can 
be sensitive to cavities with density contrast ~ 0(1) and the size (10 - 20) 
km situated not too far (< hundreds km) from the detector. The cavity can 
produce a change of the flux by about 0.1 % at most. We have found that 
exposure time during the year can reach about 10 hours. So, the detection 
of a cavity at 3a level would require a detector of about 130 Mt size. 

There are some ways to reduce the mass of detector: 

1) The exposure time can be increased if the detector moves following 
the neutrino trajectories which cross the cavity. For this, of course, some a 
priori knowledge of possible position of the cavity is needed. 

Another possibility is to perform observations in the same place during 
several years. 

2) If sterile neutrinos participate in oscillations one can use coherent ef- 
fects on scattering of neutrinos on nuclei by neutral currents in order to 
increase the cross section |2*Tj . 

Refraction effect increases with neutrino energy: e oc E. For the boron 
neutrinos we get e ~ 0.03. Moreover, the cross section of the ve— scattering 
is proportional to neutrino energy. Finally, the background conditions are 
better. So, one may think to use boron neutrinos instead of beryllium neutri- 
nos. However, the flux of boron neutrinos is 3 orders of magnitude smaller, 
and therefore total number of events turns out to be one order of magnitude 
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smaller than in the case of beryllium neutrinos. There are other problems: 
the oscillation length is now about 300 km, and therefore with boron neutri- 
nos one can probe larger structures than with Be— neutrinos Furthermore, 
due to continuous spectrum the problem appears with averaging over energies 
of neutrinos. 

Further studies are needed to show feasibility of this method to study of 
small scale structures inside the Earth and for the geological research. 
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